Capillary bridging and long-range attractive forces in a mean-field approach 
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When a mixture is confined, one of the phases can condense out. This condensate, which is oth- 
erwise metastable in the bulk, is stabilized by the presence of surfaces. In a sphere-plane geometry, 
routinely used in atomic force microscope (AFM) and surface force apparatus (SFA), it can form a 
bridge connecting the surfaces. The pressure drop in the bridge gives rise to additional long-range 
attractive forces between them. Minimizing the free energy of a binary mixture we obtain the 
force-distance curves as well as the structural phase diagram of the configuration with the bridge. 
Numerical results predict a discontinuous transition between the states with and without the bridge 
and linear force-distance curves with hysteresis. We also show that similar phenomenon can be 
observed in a number of different systems, e.g. liquid crystals and polymer mixtures. 



I. INTRODUCTION 

Capillary condensation is a classical example of a well- 
studied phenomenon in nature |]J . It is usually observed 
in a porous medium with hydrophilic surfaces, where wa- 
ter condenses out above its bulk transition temperature. 
The opposite effect, that is capillary evaporation, is ob- 
served when the surfaces are hydrophobic. 

First explained by Lord Kelvin for the case of vapor- 
water coexisting phases, capillary condensation was also 
observed in other two phase gas-liquid systems and stud- 
ied in various confining geometries 0, 0. 111 I n fact, 
capillary condensation exists in any confined system close 
to the transition point (binodal) provided that the sur- 
faces prefer one of the phases over the other and the bulk 
transition is of the first order. In this situation the con- 
densed phase is stabilized by the presence of the surfaces 
in the region of the phase diagram where the condensed 
bulk state is metastable. 

Interest in capillary condensation phenomenon has re- 
cently been renewed due to the fast development of high- 
precision force measuring devices, such as atomic force 
microscope (AFM) and surface force apparatus (SFA). 
Indeed, accurate measurements of the interactions be- 
tween surfaces separated by thin fluid films led to a dis- 
covery of unexpectedly strong and long-range attraction 
in some systems. In many cases, direct or indirect mea- 
surements proved that these attractive interactions are 
due to the formation of a capillary bridge. Below we 
mention some examples where, in our opinion, bridging 
could play an important role: 

(i) attraction between hydrophobic surfaces: Early SFA 
measurements revealed the presence of an attractive in- 
teraction between hydrophobic surfaces in aqueous so- 
lutions [a S B 13- These conclusions were confirmed 
by later AFM experiments El El El El This attrac- 
tion is much larger than could be expected for a van der 
Waals force. As a typical example, the observed range 



of force for silanated surfaces was up to 400nm with the 
maximum attraction or adhesion (normalized by radius 
of interacting surfaces) about 250mN/m |ll| ■ It has soon 
became clear that conventional theories of colloidal inter- 
actions fail to explain such a long-range effect 14, 15, 16] . 

Alternatively, it has been suggested that in some 
systems the hydrophobic attraction is associated with 
the nucleation of dissolved gas (gas-filled nanobubblcs 
or drying film) at the surfaces [17|. When the sur- 
faces become close, these form a gaseous bridge between 
them [l&UjJIifl' which generates the attractive capillary 
force. This scenario was confirmed indirectly, from the 
presence of the discontinuous steps in the force curves, 
variability in the jump separations, hysteresis between 
approach and retraction, and sensitivity to dissolved gas 
or another solut e llllll2Ll2ll . as well as from more direct 
experiments [H l^lMlIll • 

Contrary to the conclusions made on the basis of ex- 
perimental data, no consensus was achieved in the the- 
oretical description of capillary interactions. Although 
on purely thermodynamic ground cavitation is expected 
to occur between hydrophobic surfaces, both early the- 
ory [23| and later Monte Carlo study [28| predict that 
this should happen only when the separation becomes 
very small (i.e. comparable to the range of van der Waals 
force). The long-range capillary force was predicted only 
for a situation, where pre-existing hemispherical bub- 
bles (of height and radius comparable to the range of 
expected attraction) are postulated on an isolated sur- 
face [23. Contrary to this, it has been shown that the 
pre-formed nanobubbles are unstable and tend to dissolve 
very fast y(J, l3ll l32l ] . This led to speculations that the 
hydrophobic attraction is essentially a non-equilibrium 
force [13| . 

(ii) nematic liquid crystals: An exotic, but rather im- 
portant example of capillary bridgin g ha s recently been 
observed in nematic liquid crystals [331 IM . l35| . Here 
the two phases forming the bridge are both liquids with 



approximately the same density; the main difference be- 
tween them is a degree of orientational molecular order- 
ing: anisotropic molecules are orientationally disordered 
in an isotropic phase; in a ncmatic phase they are ori- 
ented along some preferred orientation, called the direc- 
tor. The fact that the two phases have different symme- 
try implies that the nematic-isotropic transition is of the 
first-order 36]. ft was confirmed 37] that, in the frame of 
the Landau-de Gennes mean-field approach, the bridge of 
the nematic phase is formed between two spherical par- 
ticles immersed in an isotropic phase, once the system is 
close to the nematic-isotropic transition. 

(iii) various mixtures: More complex in composi- 
tion, strongly confined fluids, when brought close to 
the phase separation, often have similar features in the 
force-distance profiles: the force has an unprecedentedly 
long-range attraction with a clear hysteresis on retrac- 
tion/approach. For instance, capillary bridging was ob- 
served in bicontinuous microemulsions 3{j and in a so- 
lution of polymer mixtures |3^ . 

The importance of capillary bridging as a possible driv- 
ing force in colloidal systems |J, |g, H(J, EJJ cannot be 
overestimated. The formation of the condensate between 
the particles can trigger their aggregation and, as a re- 
sult, flocculation in colloidal suspensions. Small amount 
of added fluid generates, through capillary effect, strong 
adhesive forces between the grains of granular materi- 
als mnnii. 

The direct optical measurements of the bridge param- 
eters are impossible, due to its submicroscopic size. One 
can study only the indirect consequences, e.g. force- 
distance curves in AFM or SFA experiments, and from 
there attempt to deduce some of the system properties. 
This requires theoretical modeling that fills in the gap. 

The aim of this paper is to study the general con- 
sequences of the formation/annihilation of a capillary 
bridge upon change of the thermodynamical parameters 
or geometry of the system. Considering that recent ex- 
perimental results on capillary condensation/evaporation 
are performed using AFM or SFA setups, we work in a 
sphere-plane geometry which is relevant for the above- 
mentioned techniques. We first treat the problem in 
a sharp-interface limit and derive simple analytical ex- 
pressions for the structural phase diagram and forces in 
Sec.|n] Then, in Sec. lfffl we use the Landau mean-field 
approach and solve the full nonlinear equations for the 
bridge shape numerically. In the final Section [V] using 
simple mapping of the free energies, we discuss implica- 
tion of our results to nematic liquid crystals. 



II. SHARP INTERFACE LIMIT 

To begin with, we construct a crude analytical solu- 
tion that will give us some qualitative picture of what is 
happening with the bridge when the geometry or ther- 
modynamical parameters of the system change. The ap- 
proach we shall follow was proposed by Petrov et. al. [38j. 




FIG. 1: Illustration of the studied geometry and bridge 
parametrization. Condensate of the phase (1) is stabilized 
between the apex of the sphere and the flat substrate. Axial 
symmetry is assumed throughout the paper. 



To keep the description as simple as possible, we neglect 
the internal structure of the interfaces. Then our sys- 
tem involves domains of well-defined phases separated 
by infinitely thin interfaces, so called "sharp- interface 
limit" 45]. 

A geometry typical for SFA and AFM measurements 
is depicted in Fig^ At the moment of capillary conden- 
sation the condensed phase is formed between the apex 
of the sphere and the flat substrate. For convenience, we 
denote the phase inside the bridge as a phase (I) and 
the bulk phase as a phase (2). We also assume that the 
sphere and the plane are made of the same material and 
denote it as a phase (3). Then the excess free energy, 
i.e. the change in the free energy of the system due to 
the emergence of the phase (1) (bridge) out of the bulk 
phase can be written as the sum of the bulk and surface 
terms 



AG = G b + G s 



(1) 



In the sharp-interface limit, the surface free energy G s is 
due to different surface tensions of the interfaces between 
the phases 1/3 and 2/3, plus the energy of the interface 
between the phases (1) and (2), which is proportional to 
the perimeter of the bridge 



G s = Si(ai t a - 02,3) + *SWi,2, 



(2) 



where S\ and S% are the areas of the interfaces between 
the phases 1/3 and 1/2 correspondingly. 

Additionally, we assume that the phase (1) wets the 
surfaces, i.e. the sphere and the plane are covered by 
thin wetting layers of the phase (I). This implies that the 
phase (1) has zero contact angle and therefore, according 
to the Young equation 



01,3 — 0"2,3 



-Cl,2- 



(3) 



Note that here we implicitly assumed that we are be- 
low the prewetting transition line, i.e. the layers of the 



metastable phase are thin and we can neglect their con- 
tribution to the total free energy. 

The bulk free energy is due to the difference in the 
chemical potentials of the phase-separated component in- 
side the bridge and in the bulk 



G b = VAfi, 



(4) 



where V is the volume of the bridge. 

To calculate the surface areas and the volume of the 
bridge, its circumference is approximated to a cylinder, 
instead of the correct, concave shape. Similar approx- 
imation was used before |2||. In fact, for the typical 
dimensions of the surfaces and separations of SFA ex- 
periments this approximation works very well 38] . The 
surface areas and the volume of the bridge then read 



S\ — 2nrl + irr b , 
S 2 = 2irr b (l + h), 

V = 7rr 2 (l + h)-^l(3r 2 + l 2 ), 



(5) 
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FIG. 2: Typical dependence of the excess free energy G(x,r) 
on the radial extent of the bridge for r/\ — 10. Three sep- 
arations are shown: h — 0.5A, when the bridge is stable; 
h — 0.9A, with metastable bridge; h — 1.5A, when the config- 
uration with the bridge is unstable. 



where I = r — yfr 2 — r 2 . 

It is convenient to introduce a new constant with di- 
mension of length 



A = (Ti^/A/z, 



(6) 



and scale out the variables with dimension of length by 
defining 



f = r/A, 

h = n/X, 
h = h/X, 



(7) 



We also introduce a dimensionless radial extent of the 
bridge which plays a role of the order parameter in our 
system 



= n/re [0..I]. 



(8) 



With this notations, the excess free energy JQ) can be 
written as 



AG(r, r b , h, A) = A^A 3 G(x, f, h). 



(9) 



For every fixed f, h, the minimum of G(x, f , h) gives the 
equilibrium radius of the bridge. 

In what follows it will be understood that scaling Q 
has been carried out, and we shall omit the tildes in the 
text below. 

Now that the mathematical problem is formulated, we 
first analyze the dependence of the excess free energy 
on the bridge radius r b , which is shown in Fig. It is 
similar to the Landau free energy as it occurs for phase 
transitions with a scalar order parameter x. It predicts a 
discontinuous transition between the bridged (x > 0) and 
unbridged (x — 0) configurations on approach/retraction 
of the sphere from the plane. The first-order character 



of the transition indicates possible jumps and hystere- 
sis in force-distance curves measured in AFM and SFA 
experiments [ill El HE 13 USM ■ 

Let us now localize the transition line. The equilib- 
rium radius of the bridge corresponds to the minimum 
of the excess free energy G(x, r, h). In addition, the con- 
figuration with the bridge is energetically preferable if 
G(x t ,r,h t ) < 0. Therefore, the following equations de- 
termine the transition line (binodal in a thermodynami- 
cal sense) 

dG 

— (x t ,r,h t ) = 0,G(x t ,r,h t ) = 0. (10) 

The configuration with the bridge (or without it) can 
also be metastable, when it is a local minimum of the ex- 
cess free energy. The corresponding metastability limits 
(spinodals) can be determined from the condition 



dG 



d 2 G 
r, Km) = 0, — —(x m ,r, h m ) = 0. (11) 



Now that we have a criteria for the stability and the 
metastability of the bridge, Eqns. (|10lll|) . let us have a 
look at the dependence of the bridge radius r b /r on the 
sphere-plane separation, h. If we differentiate the excess 
free energy G(x, r, K) with respect to x, we find that the 
extrema of the free energy are solutions of the implicit 
equation for the bridge radius 



\J\ - x 2 - 1 




(12) 



These solutions are shown in Fig. |3| for several values of 
the sphere radius r. As anticipated, there are two solu- 
tions for each separation h < h m (r) , which correspond to 
two extrema of the excess free energy (see Fig. |2). The 
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FIG. 3: Bridge radius rt/r as a function of the sphere-plane 
separation ft/A. The curves correspond to the different radii 
of the sphere r/X — {10, 30, 100}. Solid and dashed lines are 
binodal and spinodal correspondingly. 



larger of these two values corresponds to the minimum of 
the free energy; its value defines the radius of the stable 
(or metastable) bridge. If the separation h is larger than 
some value, h m (r), then the solution to Eq. (|12() does not 
exist. Separation h m (r) provides, therefore, metastabil- 
ity limit for the configuration with the bridge. 

Calculating the first derivative of G at x — we find 
that dG/dx(x = 0) = 2h > 0, i.e. our model predicts 
that the configuration without the bridge can be stable 
or metastable and never unstable. Both the transition 
line, Eq. (|10fl . and the metastability line, Eq. I|ll|). are 
shown in Fig. [3J 

Finally, we show the transition line (binodal) together 
with the metastability limit (spinodal) on the structural 
phase diagram, Fig. 0] It confirms that for small sphere 
radii the configuration with the bridge becomes unstable 
almost immediately, for very small sphere-plane separa- 
tions. As the radius of the sphere grows, the transition 
occurs at much larger separations, approaching h/X = 2 
for r — > oo. As a result, our model predicts that bridg- 
ing is not possible for h > 2 A. The same conclusion was 
made in Ref. |38J |. 

After the structural phase diagram is calculated, let us 
have a look at the force-distance curves. The interaction 
force can be calculated from the free energy 



T 



A^A 2 



dG(x,r, h) 
dh 



-7rr(2a; + x r), 



(13) 



and, together with equation i|12|) . provides the force- 
distance curves written in parametric form. 

This dependence is shown in Fig. together with 
metastability limits and the transition line. It quantifies 
our previous conclusion that the transition is of the first 
order: for small separations the configuration with the 
bridge is stable and gives rise to practically linear increase 
of the force with separation. At the transition point this 
configuration becomes metastable and the bridge can dis- 



FIG. 4: Structural phase diagram of the capillary bridge. For 
small r/X the bridge becomes unstable already for small sep- 
arations h/X. For large r/X the threshold separation asymp- 
totically approaches 2A. Solid (dashed) line shows binodal 
(spinodal) correspondingly. Symbols presents the results of 
the numerical minimization. Note log scale for the sphere 
radius r. 
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FIG. 5: Force-distance curves. Different symbols correspond 
to different radii of the sphere r/X = {10, 30, 100}. Solid and 
dashed lines are binodal and spinodals correspondingly. 



appear as soon as the fluctuations of the bridge radius are 
big enough to destroy it. Finally, at even larger separa- 
tions, the bridge becomes completely unstable and dis- 
appears. There is no interaction between the sphere and 
the plane from this point on and the force jumps to the 
zero value. 

On the other hand, when the sphere approaches the 
surface, the configuration without the bridge is sta- 
ble (metastable) up to the zero separation. Therefore, 
the largest hysteresis in force-distance curves on ap- 
proach/retraction is given by h m (r). 

Before presenting the numerical results, let us briefly 
recall the main conclusions. Assuming that the interface 
is infinitely sharp and the bridge has a cylindrical shape, 



we are able to see the first-order structural transition be- 
tween the configurations with and without the bridge. 
Structural phase diagram reveals that for small sphere 
radii the transition occurs almost immediately, for small 
sphere-surface separations. On the other hand, even for 
large sphere radii, bridging is not possible for h > 2<t/A/j,. 
The first-order nature of the transition leads to hysteresis 
in the the force-distance curves on approach/retraction 
of the sphere from the surface. The attractive force is 
basically proportional to the separation, with small de- 
viations close to the transition line. 



III. MEAN-FIELD APPROACH 

It is clear that even the primitive model considered 
above allows a qualitatively rich description of the tran- 
sition. However, there are still several important details 
missing, on the first place: (i) the non-cylindrical shape 
of the bridge; (ii) the finite thickness of the interface; 
(iii) the interaction between the fluid and the solid walls. 
Inclusion of these corrections makes the problem rather 
involved and does not allow anymore for a simple analyt- 
ical solution. In what follows we use numerical methods 
to tackle the problem. 



A. Free energy 

To describe the bulk phase as well as the interface 
structure, mean-field theories [4j| are often used. In 
this approach the order parameter <fi is introduced. For 
a liquid-vapour model <f> is just the density. For a bi- 
nary mixture is a composition variable, defined as 
<f> = (ni — ri2)/(ni + n 2 ), where the rii are the number 
densities of the two species. This order parameter varies 
slowly in the bulk regions and rapidly on length scales of 
the interfacial width. The unmixing thermodynamics is 
described via a free energy functional. 

In the mean-field approach the semi-grand potential of 
a binary mixture is written as |47j 



n{cj)}= fdV 



| (v^) 2 + m - A?* 



(14) 



where /(</>) is the Helmholtz free energy density of the 
mixture, while A/z is the chemical potential thermody- 
namically conjugate to the order parameter <j). 

Since the material is confined in a container in any ex- 
periment, phase separation is always affected by surface 
effects 0, 13 ■ To include them, appropriate surface terms 
responsible for the interaction of the liquid with the con- 
tainer walls are added to the free energy |4y, 1471 l4q . 
We shall note that the surface terms are very important: 
they give rise to a prewetting transition and, together 
with A/it, determine the thickness of the wetting layer. 
However, the prewetting transition occurs only in the 
vicinity of the bulk phase separation (binodal), i.e. close 



to A/i = 0, while capillary bridging can be observed in 
the whole region of the phase diagram between the spin- 
odal and the binodal. To avoid additional complications 
related to the wetting/prewetting transitions, we assume 
that our system is below the prewetting transition line 
(i.e. we have thin wetting films) and fix the surface value 
of the order parameter to its bulk value in the infinite 
system. 

The explicit form of the Helmholtz free energy f(<j>) 
varies depending on the type of mixture. However, the 
simple observation that the two phases must coexist im- 
plies that there are two minima in the free energy at the 
respective values of the order parameter. Here we adopt 
the mean-field model for a symmetric mixture |43, lH| 



a , 2 , b 4 



/W = "2^ + 4 



(15) 



Note that this is one of the simplest models to describe 
unmixing; a more realistic description would need a more 
sophisticated function, which also takes into account a 
dependence on the overall density. 

The term (V</>) is needed to provide spatial structure 
to the theory: at phase coexistence, there are two bulk 
equilibrium order parameter values <j>+ and <f>- with the 
same free energy density, f{4>+) — /(</>_). Without the 
gradient term, a structure with a large number of inter- 
faces between the </>+ and 4>- phase would be entropically 
favored. The term (fc/2) (V0) is the simplest one which 
penalizes interfaces. While this is justified near the crit- 
ical point, where interfaces are very wide and the order 
parameter varies smoothly, a more realistic description at 
strong segregation (where the interface becomes rather 
sharp) would require higher-order gradients, too. 

For the free energy (fl 51) the binodal is given by A/x = 
0. Two spinodals (metastability limits) are located at 
Afi = ±2a 3 / 2 /(27&) 1 /2. 

Finally, we should take into account the interaction 
between the solid walls and the fluid. To simplify our 
model, we have chosen an interaction that acts locally at 
the solid-fluid interface, imposing fixed boundary condi- 
tions at the walls, <j) = 4>-- 



B. Minimization procedure 

The phase-field models are very convenient because no 
explicit boundary tracking is needed. However, they are 
wasteful in terms of simulating bulk regions. One can 
remedy the situation by taking into account that the or- 
der parameter varies slowly in bulk regions and rapidly on 
length scales of the order of the correlation length, near 
the interfaces. In this situation, finite elements method 
with adaptive mesh size solves the problem of computa- 
tional efficiency: fine meshing is used only in the inter- 
facial regions; bulk regions are coarse-grained and have 
much larger size of finite elements. 

The equilibrium distribution of the order parameter <j) 
is obtained by minimizing the free energy functional i|14|) 



numerically using finite elements with adaptive meshes. 
During the minimization the square integration region 
L x L was triangulated using the BL2D subroutine |5lj . 
The function 6 is set at all vertices of the mesh and 
is linearly interpolated within each triangle. The free 
energy is then minimized using the conjugate gradients 
method [53 ] under the constraints imposed by the bound- 
ary conditions. 

A new adapted mesh is generated iteratively from the 
previous minimization. The new local triangle sizes are 
calculated from the variations of the free energy, in order 
to guarantee a constant numerical weight for each mini- 
mization variable |53j ]. The final meshes, with a minimal 
length of ~ 10~ 3 , had around 10 4 minimization variables. 

In order to obtain both stable and metastable con- 
figurations, we used different types of initial conditions 
including configurations with the capillary bridge and 
without it. To select stable solutions, we calculated the 
grand potential of the mixture, O, for both stable and 
metastable solutions and chose the solution with the low- 
est grand potential. This allows the accurate determina- 
tion of the phase diagram. 



IV. NUMERICAL RESULTS 
A. Free interface 

To test the minimization algorithm and to tune the 
minimal mesh size we first consider a system without a 
spherical particle and a surface. We impose the boundary 
conditions in such a way that a flat interface sets in the 
system, located at z = 0, perpendicular to the z axis. 

In this geometry the order parameter 6 depends only 
on the z coordinate and the interface profile is a solution 
to the following Euler-Lagrange equation 



dz 2 06 



(16) 



together with the boundary conditions </>(±oo) = ±<fo, 
where ±<^>& are the homogeneous bulk solutions, given by 
df/86 = 



6± = ±6 b = ±y/a/b. (17) 

The solution to the boundary problem (|16fl yields 



6 — &b tanh (z/t;) . 



(18) 



where £ = y / '2k/a is the correlation length, which serves 
as a measure of the thickness of the interface. 

With the solution (|18fl we go back into the free en- 
ergy 1)14(1 . integrate over the z coordinate and obtain the 
interfacial tension a, defined as the excess free energy of 
the interface 



= k 



dz 



-Mlz 



(19) 
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FIG. 6: Order parameter profile of a free interface. Cir- 
cles: results of the numerical minimization; solid line: fit to 
Eq. 1181 with £ = 0.02235. Inset illustrates the adaptive 
mesh: smaller triangles are used in the interface region with 
strong variation of the order parameter. 



For our numerical minimization we used the following 
parameters of the potential (|14|) : a = b = 4, k = 0.001. 
They yield the correlation length £ » 0.02236 and the 
interfacial tension a « 0.06. 

The order parameter profile for this set of parameters 
is shown in Fig. [fj] As expected, the order parameter 
changes between two bulk values ±0h = ±1 and the 
width of the interface is of the order of the correlation 
length £. A fit of the results of the numerical minimiza- 
tion to Eq. I|18|) yields £ = 0.02235, practically indistin- 
guishable from the exact value y/2kja. 

The interface thickness £ sets the smallest length-scale 
in our system and, therefore, defines the minimal size 
of the mesh. We found that a minimal length of the 
final meshes £ = 10~ 3 , which is about 20 mesh points 
per interface, is accurate enough to recover the interface 
structure. 

The correlation length £ and the interfacial tension a 
are two parameters of our system which can be experi- 
mentally measured. In addition, the interfacial tension 
enters the sharp-interface description presented above. It 
is, therefore, convenient to rewrite the excess free energy 
in terms of these two parameters. This shall help us to 
compare the numerical results to the sharp-interface limit 
as well as to the existing experiments. 

Performing the substitution ip — &/ 4>b we obtain the 
semi-grand potential density in the form 



3a 

U) = — 
8£ 



-if 



2^ 2 + ^ 






As before, we introduced a constant 



A 



(20) 



(21) 



2A/l^ 6 

with dimension of length. Note that the length A has 
(approximately) the same value as in the sharp-interface 



limit. To prove this, consider the situation when £ — > 0. 
Then the interface is sharp and the part of the excess free 
energy which scales as the bridge volume can be written 
as 



f ex = [A//(^+ - <M + /(<M - ftf+)]v, 



(22) 



where V is the volume of the bridge. For small chem- 
ical potential differences <f>± ~ ±</>6i where cf>b is the 
value of the order parameter when A/j, = 0. Then f ex 
can be rewritten as 2Afi(j)bV. Comparing f ex to the ex- 
pression for the bulk free energy in the sharp interface 
limit, Eq. (0J, we obtain the definition of A in the form 
of Eq. (EIJ. 

Now that the semi-grand potential density is written in 
a dimensionless form it is clear that there is an additional 
length-scale in the problem, when it is formulated in the 
phase-field approach. This length-scale is the correlation 
length £ or, alternatively, the width of the interface be- 
tween the two phases. 

Three more variables with dimension of length are 
present in the problem. Two of them, sphere radius r and 
sphere-plane separation h, fix the geometry of the system. 
For both SFA and AFM experiments r > h. The third 
variable, A, specifies the thermodynamic state of the sys- 
tem. The thickness of the interface £ validates the sharp- 
interface limit: it is applicable only when £ -C A, h, r. 




FIG. 7: Cross-section of a capillary bridge for several surface- 
tip separations: ft/£ = {0,7.16,11.63,13.42}. Sphere radius 
r/£ = 22.36, A/£ = 13.42. b) illustrates the final mesh; c) 
shows the interface contour lines, which are parallel to the 
surface of the bridge. 



B. Sphere-plane geometry 

We have thus seen that the interface thickness £ sets 
the length-scale in our system. Therefore, three dimen- 
sionless ratios uniquely specify the state of the system: 
r/£, h/£, and A/£. In what follows we fix the thickness of 
the interface to £ = 0.02235; the corresponding constants 
of the potential are discussed in Section IIV Al The rest 
of the parameters with dimension of length (i.e. r, h, A) 
are given in units of £. 

Let us first look at the situation when the thermody- 
namic parameters of the system (i.e. A, in addition to 
£) are fixed, but the geometry is changing. For these 
calculations we used the chemical potential difference 
A[i = 0.1, which yields A/£ « 13.4. 

The functional 114|) was minimized for a range of 
sphere radii, r/£ 6 [0.5..200] with the boundary condi- 
tions 4> = </>- a t the sphere and the plane surfaces. Note 
that the smallest size of the sphere is naturally set by 
the phenomenological description we use: it is not ap- 
plicable on length-scales much smaller than the interface 
width £. The upper value is limited by the computa- 
tional capacities: further increase in the sphere radius led 
to significant slowdown of already rather time-consuming 
calculations. 

Typical cross-sections of the bridge are shown in Fig.0 
For small sphere-plane separations the wetting film, cov- 
ering the sphere and the plane, becomes thicker next to 
the sphere, Fig. [7f a). On retraction of the sphere from 



the plane the bridge is formed and then gradually thin- 
ners, see Figs. EJb) and (c). Finally, at even larger sep- 
arations, the bridge becomes unstable and disappears. 
The configuration with the sphere and the plane, wetted 
by a thin wetting layer, Fig.[7{d), becomes energetically 
favorable. At these distances the sphere and the plane 
no longer interact with each other. 

To investigate the nature of the structural transitions 
between the bridged and unbridged configurations, we 
plot the excess free energy as a function of the sphere- 
plane separation h in Fig. |H1 

The free energy profiles quantify the conclusions made 
above: when the sphere is far from the plane, the con- 
figuration without the bridge is the only stable one. The 
excess free energy does not change with separation, i.e. 
there is no force acting on the sphere. When the sphere 
moves closer to the plane, this configuration becomes first 
metastable, and then unstable. The bridge is formed and 
stays until the sphere touches the surface. 

On the other hand, on retraction, the configuration 
with the bridge becomes first metastable and then un- 
stable. Jump in the free energy profile points on dis- 
continuous (first-order) structural transition between the 
configurations with and without the bridge. Due to 
the first-order nature of the transition, the configura- 
tion with the bridge can be kinetically stabilized and 
the bridge can, in principle, disappear at any position 
where it is metastable. In fact, this explains the vari- 
ability in the jump-in/out distances observed in experi- 
ment [HEiil. 
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FIG. 8: Excess free energy of the configuration with the 
bridge vs sphere-plane separation. Different curves corre- 
spond to different sphere radii. Inset illustrates change in 
the radial extent and shape of the bridge when the sphere 
radius changes. The inset configurations are taken at a fixed 
sphere-surface separation, h/£ — 8.94. 



On the face of it, the scenario predicted by the phase- 
field approach is similar to the previous, sharp-interface 
limit one. However, in addition to quantitative correc- 
tions, it also predicts qualitatively different results for 
small r/£. Indeed, as it is seen from the free energy 
profile for r/£ = 0.9, there is no jump in the excess free 
energy. Instead, we have a smooth cross-over between 
the bridged configuration and configuration without the 
bridge. This points to a critical point r c below which the 
bridge will no longer disappear discontinuously The dif- 
ference between the phase-field and the sharp-interface 
approaches for r/£ w 1 is, of course, expected, since the 
infinitely sharp interface assumption breaks down when 
the radius of the sphere is of the order of the correlation 
length £. It is then tempting to explain the absence of 
a jump in computer simulations [28j . where the system 
sizes are typically limited to thousands of angstroms. 

To calculate the transition line, we fitted the free 
energy profiles with a simple polynomial dependence 
(quadratic polynomial turned out to be a reasonable ap- 
proximation) and found where it intersects the Q = 
axis. This allows accurate determination of the tran- 
sition line. The results are shown in Fig. 0] together 
with the structural phase diagram obtained in the sharp- 
interface limit. The agreement above the critical point 
is reasonable, taking into account that we made rather 
crude approximations about the shape of the bridge in 
the sharp interface limit. 

We have thus seen that the excess free energy profiles, 
scaled out by the sphere radius, are accurately fitted by 
a quadratic polynomial. For large r/£ these profiles have 
convex shape pointing that the attractive force, which is 
the first derivative of the excess free energy with respect 
to the sphere-plane separation, decreases linearly with 
the increase of the sphere-plane separation. It can be also 



FIG. 9: Force-distance curves calculated for several dimen- 
sionless sphere radii r/£. Inset shows the slope of the curves 
a vs the sphere radius r. 



seen that the free energy profiles, when scaled out by r, 
are almost parallel to each other, especially for r/£ — > oo. 
This means that the slope of the force, scaled by the 
sphere radius r, approaches asymptotically a constant 
value when r/£ — > oo. For smaller r/£ the force profiles 
become first linear in h and then start to curve inwards, 
changing their shape to concave. As a result, the force is 
independent of h for some particular sphere radius; for 
even smaller r/£ it increases with the increase of h. 

This is illustrated in Fig. where we show the force- 
distance curves calculated for A/£ = 13.4. The inset 
shows the slope a obtained from the linear fit 



T = Tq + a^arh. 



(23) 



Note that here both r and h arc dimcnsionless, i. e. given 
in units of £. The fit confirms that for r/£ > 10 we have 
a constant slope a ~ 0.5, i.e. the force simply scales 
with the sphere radius r, like in the Derjaguin approxi- 
mation J54J. 

We studied two more situations, with A/£ = 134 ^> 1 
and A/£ = 1.34 ~ 1. The transition lines for both cases 
are also shown on the structural phase diagram, Fig. 0] 
For A/£ = 134 the agreement with the sharp interface 
limit is remarkably better. At the same time, it is clearly 
worse for A/£ ~ 1: the transition line is shifted toward 
bigger separations and ends in a critical point for rather 
large r/X. There are two reasons for this. First, the shape 
of the bridge is, in fact, a surface of revolution whose 
mean curvature is constant. These surfaces are known as 
Delaunay surfaces, or onduloids 55]. A defines the "ra- 
dius of curvature" of the bridge cross-section through its 
symmetry axis. If A 3> r the bridge is well approximated 
by a cylinder. This explains good agreement with the 
sharp interface limit for A/£ = 134. On the other hand, 
when A ~ £ the radius of curvature is of the order of the 
interface thickness, i.e. the sharp interface limit cannot 
be used, even if r, h 3> £. 



In spite of the complications in the structural phase 
diagram for A ~ £ and r ~ £, the main conclusions of our 
numerical studies remain: the transition between bridged 
and unbridged configurations is of the first order; the 
transition line ends in a critical point located at r c /£ ~ 1. 
For large spheres, r/£ > 10, the force is a linear function 
of the sphere-plane separation h with a slope proportional 
to the sphere radius r. 



V. EXAMPLES AND DISCUSSION 

In the previous sections a theoretical and numerical 
analysis of the bridge formation and stability has been 
developed for a simple binary mixture. Here we discuss 
implication of our results for several different systems. 

We begin with nematic liquid crystals, mentioned in 
Introduction. From the phenomenological point of view, 
the system is described by the Landau-de Gennes free 
energy |5q 



F{Q}= J(f b + fe)dV + J f s dS 



(24) 



where /& is the bulk free energy density, f e is the elas- 
tic free energy density and f s is the surface free en- 
ergy. Within a mesoscopic approach the minimum of 
the Landau-de Gennes functional F{Q} gives the equi- 
librium value of the tensor order parameter Q. 

Symmetry arg uments yield for the local bulk free en- 
ergy density J5|| |53 



f b = aTrQ 2 - bTiQ 3 



[TtQ 2 



(25) 



where a is assumed to depend linearly on the tempera- 
ture, while the positive constants b, c are taken tempera- 
ture independent. 

It is convenient to scale out the variables by defining 



Qij = 6c/bQ i: j, 
h = 24 2 c 3 /6 4 /b- 



(26) 



It will be understood that such scaling has been carried 
out, and we shall omit the overbars in the text below. 

We also introduce a dimensionless temperature r by 
defining 



t& 2 /24c. 



(27) 



The elastic free energy density can be written as [57 



1 dQjj dQjj 

2 dx k dx k 



2 dxj dxk 



where the constants L\ and L 2 are related to Frank- 
Oseen elastic constants by K\\ — K33 = 9Q 2 (Li + 
L 2 /2)/2 and if 22 = 9Q 2 Li/2 and Qb is the bulk nematic 
order parameter. The sign of L 2 defines the preferred 
orientation of the director at the NI interface. L 2 > 
(Z/2 < 0) favors planar (perpendicular) anchoring 58] . 



Now that the full free energy is specified, we perform 
several simplifications. To begin with we neglect the non- 
uniformity of the director distribution inside the bridge. 
This assumption is quite reasonable for the studied geom- 
etry: we assume homeotropic (perpendicular) anchoring 
of the director at the surfaces and planar (parallel) an- 
choring at the nematic-isotropic interface. Second, we 
assume that the nematic phase is uniaxial and neglect 
the biaxiality of the interface. The full problem has been 
considered recently by Stark et. al. 37]. 

For a uniform uniaxial nematic (Qu = Q, Q22 — 
-1/2Q) the free energy l|25(l takes the form 



133 



f b = tQ 2 ~ 2Q 3 + Q 4 



(29) 



Free energy (|29|l predicts that the nematic state is stable 
when t < tni = 1 with a degree of orientational order 
given by 



Q b = - A 1 



1- -T 

9 



(30) 



The same expression provides us with the metastability 
limit of the nematic phase on heating at r* = 9/8. 

The solution which corresponds to a free interface is 
then given by the de Gennes Ansatz [58J . 



1 



1 + tanh ■ 



£ 



(31) 



where £ is the nematic correlation length when the direc- 
tor is parallel to the interface 



e 



8Lic 



L 2 



(32) 



A linear transformation of the order parameter, Q = 
(1 + 4>)/2, reduces the Ansatz <|31ll to the familiar 
form l|18Jl . The same substitution transforms the free en- 
ergy of a uniaxial nematic to the free energy of a binary 
mixture 



h 



2(3 - 2r)0 2 + 8(r - 1)0]/16. (33) 



Close to the nematic-isotropic transition, tni = 1, the 
r-dependence of the coefficient at </> 2 is not important, 
and the free energy has exactly the same form as H21(l 
provided that 



A = 



£ 



6(r-l) 



(34) 



This is an interesting result: recalling that the nematic- 
isotropic transition occurs at tni = 1 and the superheat- 
ing temperature of the nematic phase r* = 9/8 we see 
that A > 4/3£, i.e. A spans the whole range of values, 
from the interface thickness £ at the superheating tem- 
perature, to the system sizes, set by the sphere radius r 
and the sphere-plane separation h, and diverges at the 



10 



nematic-isotropic transition. This implies that, depend- 
ing on the temperature, one can have either discontinu- 
ous (first-order) or continuous transition between bridged 
and unbridged configurations. Moreover, by changing the 
temperature of the liquid crystal it is possible to control 
the range and the strength of the attractive force between 
the particles in colloidal suspensions. 

Another 'liquid crystalline' example for which our re- 
sults are relevant is drag of colloidal particles by the 
nematic-isotropic interface. Recent experiments | 59j re- 
vealed that the nematic-isotropic interface is able to drag 
colloidal particles and, as a result, form spatial patterns 
in the liquid crystal cell. Theoretical reexamination of 
the experimental data J23 demonstrated that the force 
exerted on a colloidal particle is proportional to the pen- 
etration depth, and the slope of the force scales as the 
particle size r. This points to rather unexpected map- 
ping: here the bridge connects the colloidal particle and 
the interface itself, since the position of the latter is fixed 
by a small temperature gradient present in the system. 
The straightforward implication of our results, which is 
also confirmed experimentally |59| , is that it is more dif- 
ficult for the interface to capture small particles. Indeed, 
the attractive force from the side of the interface scales 
with radius. In addition, the barrier becomes smaller, 
completely vanishing at the critical radius r c . 

Finally, we would like to mention that the standard 
Flory-Huggins expression for a free energy of an incom- 
pressible polymer mixture [6ll l6^ | can be expanded 
in powers of the polymer volume fraction c\>. Close to 
the critical point it has the same form as l|15|) with 
a = -k B T/N+ X /2, b = k B T/{3N), where x is the Flory- 
Huggins parameter, N is the effective chain length. This 
means that the whole picture of bridging, described here, 
can be adapted for the case of incompressible polymer 
mixtures. In fact, the capillary-induced phase separa- 
tion in binary polymer solutions has been experimentally 
observed [39j and theoretically studied using mean-field 
lattice modeling |63| . 

In many ways the interface acts like a membrane un- 



der tension. In our particular situation the formation 
of a bridge between the surface and a particle resembles 
the process of wrapping a colloidal particle by a mem- 
brane |64ll65| . A detailed comparison requires, however, 
further analysis. 

As a final remark, we remind that, by using rigid 
boundary conditions, we will not observe possible prewet- 
ting transitions which have been studied in detail for 
simpler geometries 0, IS |66|. These transitions nor- 
maly occur close to the bulk phase separation, i.e. when 
Afi w 0, or A > 1. The corresponding part of the phase 
diagram (r/X,h/X -c 1) can, therefore, be affected by 
these effects. On the other hand, when r — > oo, we re- 
cover the capillary condensation in a slab geometry |2|. 
In fact, the critical point has similar to the capillary con- 
densation origin: it occurs when the interfacial width is 
comparable to the sphere-plane separation, i. e. to the 
width of the capillary. 

VI. CONCLUSIONS 

In conclusion, we have presented a detailed study of 
interactions between the sphere and the plane due to for- 
mation of the capillary bridge between them. We have 
shown that the structural transition between the configu- 
rations with and without the bridge is discontinuous and 
ends in a critical point for small particle radii. We have 
also demonstrated that the attractive force is long-ranged 
and is basically proportional to the sphere-plane separa- 
tion and possible hysteresis on approach/separation. 
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